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Abstract 



We use covariant techniques to describe the properties of the Godel universe and then 
consider its linear response to a variety of perturbations. Against matter aggregations, we 
, find that the stability of the Godel model depends primarily upon the presence of gradients 

^ ' in the centrifugal energy, and secondarily on the equation of state of the fluid. The latter 

, dictates the behaviour of the model when dealing with homogeneous perturbations. The vor- 

ticity of the perturbed Godel model is found to evolve as in almost-FRW spacetimes, with 
some additional directional effects due to shape distortions. We also consider gravitational- 
• wave perturbations by investigating the evolution of the magnetic Weyl component. This 

, tensor obeys a simple plane-wave equation, which argues for the neutral stability of the Godel 

' model against linear gravity-wave distortions. The implications of the background rotation 

Q I for scalar-field Godel cosmologies are also discussed. 

^■ 

^ ; PACS numbers: 04.20.Cv, 98.80.Jk 

• r-j . 1 Introduction 
X 

. The Godel universe is an exact solution of the Einstein held equations, which is both stationary 

and spatiahy homogeneous [1]. The model, which is of Petrov type D, is also rotationally sym- 
metric about each point and contains a perfect-fluid matter source whose 4-velocity is a Killing 
vector [2]. Godel's universe is known for its unusual global properties. The most intriguing 
among them is the existence of closed timelike curves, which violate global causality and makes 
time-travel theoretically possible in this space-time. Hermann Weyl had flrst suggested in 1921 
that time travel might occur in general relativity [3] as a consequence of "very considerable 
fluctuations in the space-time metric" but he believed that the fluctuations "necessary to pro- 
duce this effect do not occur in the region of the world in which we live". It was clear that 
this causal anomaly could arise for some distorted geometry and hence for some distribution of 
mass and energy. Kurt Godel showed that the required distribution could be extremely simple 
and his discovery initiated the study of the global properties and causal structure of Einstein's 
equations [4]. At flrst, attention was almost entirely focussed upon the properties of the equa- 
tions of motion in the Godel universe [5]- [7] and an interesting overview of these investigations 
has recently been given by Ozsvath and Schiicking [8]. Subsequently, Godel's solution has also 
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triggered a considerable amount of work on rotating solutions of the Einstein field equations 
and rotational cosmological models (sec [9]- [13] for a representative list). For further details and 
an extensive discussion the reader is referred to the recent review articles by Krasihski [14]-[16] 
and Obukhov [17]. Although the Godel spacetime is not a realistic model for our universe, it is 
an important theoretical laboratory for investigating a range of global properties of the space- 
time structure in different gravity theories. Recent studies of quantum computation [18, 19] have 
shown that the presence of closed timclikc curves in spacetime provides for a new physical model 
for quantum computation in compact regions. A quantum computer with access to closed time- 
like curves can solve NP-complete problems with only a polynomial number of quantum logic 
gates. A scries of studies of the possibility of closed timelike paths in string theories [20]- [26] 
have also drawn on the insights gained from the detailed study of the Godel spacetime and its 
generalisations to arbitrary space dimensions [27, 28], considering the constraints of possible 
holographic principles [29], and investigating the accessibility of information in different parts 
of spacetime [30]. 

The Godel universe is a member of the family of homogeneous spacetimes together with 
the Einstein static and de Sitter universes. In recent years close attention has been paid to the 
stability of the de Sitter universe and its physical significance for the consequences of inflation in 
the early stages of our universe. Under physically realistic conditions on the energy-momentum 
tensor of matter, all perturbations to isotropy and homogeneity will be seen to fall off expo- 
nentially rapidly within the event horizon of a geodesically moving observer in the de Sitter 
universe [31]- [33]. This result, in its different technical expressions, is known as the cosmic no 
hair theorem. The stability of the Einstein static universe has also been examined in a covari- 
ant fashion and reveals a more intricate dependence on material content than is conventionally 
reported [34, 35]. We will complement these detailed studies of the stability of the de Sitter 
and Einstein static homogeneous spacetimes with a covariant study of the stability of the Godel 
universe. Despite the past interest in the causal structure of the Godel model there have been 
few perturbative studies of its properties and these have been confined to the investigation of 
the effects of the global rotation on the evolution of scalar density perturbations in the papers 
of Silk [36, 37]. Here, we investigate the stability with respect to scalar, vector and tensor per- 
turbation modes using the gauge covariant formalism of Ellis and Bruni [38] . In the process we 
show that the background vortical energy contributes to the gravitational pull of the matter, 
while it gradients add to the pressure support. The balance between these two agents effectively 
determines the stability of the Godel universe against matter aggregations. By examining the 
rotational behaviour of the perturbed model, we explain how shape distortions can increase or 
decrease its overall rotation. We also introduce a set of linear covariant constraints, which isolate 
the pure tensor perturbations, and analyse their propagation on the rotating Godel background. 
Our results argue for the neutral stability of these gravitational-wave distortions. Finally, we 
consider some of the effects of spatially homogeneous perturbations and discuss the compatibility 
of scalar fields with the symmetries of the Godel spacetime. 



2 Covariant characterisation of the Godel universe 
2.1 The irreducible kinematical variables 

The covariant description of the Godel universe, with respect to its irreducible kinematical 
quantities, has been given in [39] (see also [4, 40]). Here, we will briefly introduce and extent 
this description and also present the associated constraints. 
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Relative to a timelike 4- velocity field Ua (normalised so that UaU°- = — 1) that is tangent to 
the worldlines of the fundamental observers, the Godel spacetime is covariantly described by [39] 

Q = = ita = CTab , and Wa 7^ . (1) 

Therefore, with the exception of the vorticity {uja), the rest of the kinematical variables, namely 
the volume expansion (6), the shear (aab) and the acceleration vanish identically. We note 
that Vba;a = 0, so ensuring that the vorticity vector associated with the 4-velocity field Ua is 
covariantly constant. 

2.2 The twice-contracted Bianchi identities 

The stationary nature of the Godcl solution means that there are no propagation equations: they 
have all been transformed into constraints. For example, on using the twice-contracted Bianchi 
identities we obtain the standard conservation laws for the energy and momentum densities. 
When applied to the Godel spacetime, the latter yield the constraints 

fi = and BaP = 0, (2) 

where /i and p are the matter density and isotropic pressure respectively.^ Here, an overdot 
indicates differentiation along (e.g. fi = u"''Val^)- Also, = hJ'Vb is the covariant derivative 
operator orthogonal to Ua, and hab = gab + UaUb is the associated projection tensor. 

2.3 The Ricci identities 

Covariantly, the kinematic evolution is determined by a set of three propagation equations and 
three constraints, all of which derive from the Ricci identities V[aV;,]Uc = Rdcbau'^i where Rabcd is 
the spacetime Riemann tensor. In Godel's universe, the Raychaudhuri equation, which describes 
the volume evolution of a fluid element, reduces to 

iK(/i + 3p) - 2a;2 - A = , (3) 

with K = SttG and uj'^ = uJaUJ^- Note how the rotation balances the gravitational attraction of 
the matter, as well as that of the (negative - see Sec. 2.5) cosmological constant. Thus, in the 
Godel model the vorticity has assumed the role played by the (positive) cosmological constant 
in the Einstein static universe. Of the two remaining propagation equations, the shear evolution 
formula takes the form 

Eab + i^ia^b) = , (4) 

where Eah is the electric component of the Wcyl tensor.^ The vorticity propagation equation, 
on the other hand, is trivially satisfied. Two of the three kinematical constraints, the shear- 
divergence and the vorticity-divergence, are also trivially satisfied, while the gravito-magnetic 
constraint leads to 

Hab = Q, (5) 

^Originally, Godel's solution was given for dust (i.e. /i 7^ 0, p = 0, A 7^ 0). However, through the transformation 
IJ, ^ n' = + p and A A' = A + up, the Godel spacetime can be reinterpreted as a perfect-fluid model. 

^Angled brackets are used to indicate orthogonally projected vectors and the projected, symmetric and trace- 
free part of second rank tensors. 
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with Hab being the magnetic counterpart of It should be emphasised that the latter is not 
covariantly constant (i.e. V c^ab 7^ 0)- Instead, one can use Eq. (4) to show that 

Eab = Q = T},Eab , (6) 

in agreement with the stationary nature and spatial homogeneity of the Godel spacetime. 
Clearly, result (6b) also guarantees that D^'Eab = = curlE^a^. 

2.4 The Bianchi identities 

The Bianchi identities provide two pairs of propagation and constraint equations for the con- 
formal curvature, which is monitored through the electric and the magnetic parts of the Weyl 
tensor. When applied to the Godel spacetime, the £^-equation gives 

E^/eb)cdU^' = 0, (7) 

while the if-equation is trivially satisfied. Note that, when (4) is taken into account, the above 
constraint also becomes trivial. On the other hand, the divEab and diviJ^b constraints associated 
with the Bianchi identities lead to 

Bafi = (8) 

and 

+ p)u;a + 3EabU^ = , (9) 

respectively. Results (2a) and (8) combine to ensure that = 0, while from Eq. (3) we obtain 
p = 0. In other words, both the energy density and the isotropic pressure of the matter that fill 
Godel's universe are covariantly constant quantities. 

2.5 Further constraints 

Additional constraints are obtained by contracting (4) along Ua and substituting the result into 
Eq. (9). Then one finds that 

K(^+p) = 2L<J^ (10) 

which also measures the total inertial mass of the Godel universe. On using this result, we can 
recast Eq. (3) as 

K{f^-p) = -2A, (11) 

to guarantee that A < as long as p < /x. The cosmological constant vanishes only when 
the fluid has a maximally stiff equation of state, with p = fx. Thus, for conventional matter 
sources, the Godel spacetime has non-positive A. This restriction is relaxed when dealing with 
the Newtonian analogue of the Godel universe (see below). From (10) and (11) it becomes clear 
that, given the equation of state of the matter, only one of /x, a; or A is needed to determine the 
other two. Finally, constraints (10) and (11) combine to give 

+ A-oj'^ = 0, (12) 

which is the Godel analogue of the Friedmann equation. 

The generic rotation of the Godel universe means that the fluid flow lines are not hypersurface 
orthogonal and, therefore, there are no integrable spatial sections. Nevertheless, one can still 
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employ the generalised Gauss-Codacci equation to evaluate the orthogonally projected Ricci 
tensor. Applied to the Godel model, the latter reads 



Uah = "^{t^H + - U?)hab , (13) 

which by means of constraint (12) ensures that TZab vanishes. 

2.6 The Newtonian analogue 

There is a simple Newtonian counterpart to the rotating Godel universe in the case of zero pres- 
sure that has been explored by Ozsvath and Schiicking [41] (see also [42] ) following Godel [43] . 
If we consider a pressureless fluid of constant density ^ rotating rigidly with constant angular 
velocity ui about the z— axis so v = {—ujy,iox,0) then, by integrating the Euler equation 

Ov 

— + v-Vv = -V^, (14) 
we find that the Newtonian gravitational potential is 

$ = X(^' + y')- (15) 

where now V indicates ordinary partial derivatives. This solves the Poisson equation with 
cosmological constant term (c = 1), 

V^^ + A = ^Kii, (16) 

if 

2cj^ = iK/i-A. (17) 

Wc sec that this Newtonian result coincides with the relation (3) for the Godel solution in the 
p = case and, in contrast to general relativity, can be satisfied even for A = and some A > 
values. In general relativity the additional geometrical constraint (11) leads to the stronger 
relation, 

= i/t/x , (18) 
for Godel's solution, which therefore requires A < 0. 

2.7 The closed timelike curves 

The most intriguing property of the Godel spacetime is the existence of closed timelike curves, 
which makes time travel a theoretical possibility and has ensured an enduring interest in the 
model [44]- [47]. To demonstrate the presence of such closed timelike worldlines, consider the 
line element of the original Godel solution (i.e. with p = and metric signature -2) written in 
cylindrical coordinates (r, z, (p) 



ds^ = 4a2 



- dr^ - dz^ + (sinh-^ r - sinh^ r) d^^ + 2^2 sinh^ rd^dtl , (19) 



with 1/a^ = K/J, = — 2A [1]. Then, introducing the coordinate transformation t ^ i = 2at, 
r ^ f = 2ar, z ^ z = 2az and cp ^ (p = (p, transform the above line element into 



ds^ = dt'^ - dr^ - dz^ + 4a^ 



sinh^ (|^) - sinh^ (^^) d<p^ + 4V2a sinh^ (|^) d^dt , (20) 
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where the tildas have been suppressed. Clearly, if for some r = ro the quantity sinh^(r/2a) — 
sinh^(r/2a) is positive, the circle defined hy r = rQ, t = = z will be timclikc everywhere [1]. 
Thus, on using (10) with p = 0, the condition r > 21n(l + \/2)/y^K/I for the existence of such 
closed timelike curves reads 

,>«,,^§Mi±i^, ,21) 

where Rq denotes the radius of the observer's "causal region". Hence, Rq ^ oo as ^ 0, which 
means that the weaker the rotation of the model the more "remote" the closed timelike curves 
become. Alternatively, one might say that the faster the Godel universe rotates the smaller its 
causal region becomes. Note that one can include the fluid pressure in Eq. (21) by simply using 
the transformation /j, ^ p, = fx + p. 



3 The perturbed Godel universe 

3.1 Conservation laws 

For a perfect fluid, the energy density and the momentum density conservation laws are given 
by the nonlinear expressions 

A = -(m + p)0, (22) 

and 

in + p)Ua = -BaP , (23) 

respectively. These formulae also hold after we have linearised about the Godel background. 
In that case, however, we can substitute the zero-order relation p = 2a;^ on the right-hand 
side of (22). Also, when dealing with a barotropic perfect fluid (i.e. for p = p{lJ.)), we have 
DaP = CgDa/x, where Cg = dp/d/x is the square of the adiabatic sound speed. 

3.2 Density perturbations 

Consider a general spacetime filled with a single barotropic perfect fluid. The nonlinear evolution 
of density inhomogeneities is monitored by the system [38] 

= W@Va-{aab-UJab)V'' -{l + w)Za, (24) 

-^(a) = -lQZa-{aab-u;ab)2^ -^KiJ,Va + a^Ua + aAa-2aDa{a'^ -Lo"^), (25) 

where Va = {a/fijDaH, Za = aDa& and w = p/fihj definition.^ Also 

?fi = Kfi-le^-2{a'^-Lo'^) + A + A, (26) 

defines and Aa = D,jA with A = V^iia- The 4-acceleration satisfies the momentum-density 
conservation law, which for adiabatic perturbations is given by Eq. (23). The projected vectors 
T^a and Za respectively describe local inhomogeneities in the matter energy density and in the 
volume expansion [38]. Note that in an exact Godel universe Da/U = = D(j0. Hence, V^, and Z^ 
vanish identically to zero order and therefore both satisfy the gauge-invariant requirements [48] . 

^The background scale a, which has been introduced to make V and Z dimensionlcss and appears in Eqs. (24), 
(25), is a characteristic dimension of the Godel universe corresponding to the radius of the smallest closed timelike 
curves [1, 37]. 
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The scalar ^ plays an important and subtle role via its coupling with the 4-acceleration in 
Eq. (25). In particular, the sign of 3? determines the effect of the pressure gradients associated 
with Ua- When 5ft is positive these gradients will add to the gravitational pull of the density 
inhomogeneities. Otherwise, they will contribute to the pressure support. The subtlety in 3fi is 
that the vorticity adds to the effect of the ordinary matter, whereas the shear opposes it. This is 
a purely relativistic effect, with no known Newtonian analogue. Note that in exact FRW models 
3f? = 0, which explains why the aforementioned behaviour has passed relatively unnoticed [38]. 

Linearising Eqs. (24), (25) about the Godel universe and keeping only the zero-order com- 
ponent of (26) we obtain 

■Da = CVab'D'' - (1 + w)Za , (27) 

Za = COabZ^ - i«M^a + + aAa + 2aBaCo'^ , (28) 

with 

K = + 2cj2 + A = 3a;2 > , (29) 

since k/j, — lo^ + A = in the exact Godel spacetime. Also, A = D^iia to first order and Ua is 
still given by Eq. (23). Then, the system (27), (28) reads 

Va = COab'D'' - {l+w)Za, (30) 

Za = LOab^^ -ki^fiVa- - D„p - ° D^(D^p) + 2aD^a;^ (31) 

IJ,{1+W) ii{l+w) 

where oj^ = /«/x(l -|- w)/2 to zero order (see Eq. (10)). Compared to the perturbed FRW case, 
the nonzero rotation of the Godel background has added four extra terms to the system (30), 
(31). They are the first term on the right-hand side of (30), and the first, third and the last 
term on the right-hand side of (31). Note the second of the three vorticity terms in Eq. (31), 
which is proportional to and adds to the overall gravitational pull. This relativistic effect 
seems to suggest that rotational energy also has "weight" . Indirectly, it also seems to favour the 
de Felice and the Barrabes et al interpretation [49, 50] of the Abramowicz-Lasota "centrifugal 
force reversal" effect [51, 52]. The third vorticity term in (31), on the other hand, is triggered by 
inhomogeneities in and will also appear in a Newtonian treatment. It can resist the collapse, 
thus mimicking the effects of the ordinary pressure gradients. As we shall see below, the stability 
of the Godel universe depends crucially on which of these two terms dominates over the other. 

Density perturbations are related to gradients in the vorticity and the curvature by means 
of the Gauss-Codacci equation, according to which the linearised projected Ricci scalar is given 
by TZ = 2(«;/x — uj^ + A). The above immediately implies the linear relation 

K^iVa = aBaiv^ + iaD„7^ , (32) 

between Dato^ and DaTZ- 



3.3 Two alternative types of perturbations 

Following (32) we may consider two types of density perturbation, first by setting BaTZ = 
and then by assuming that Dao;^ = 0. For simplicity, we will label the former isocurvature and 
the latter perturbations under rigid rotation. In general one does expect these two perturbation 
types to propagate independently. However, this approximation scheme simplifies the system 
(30), (31) and we can investigate the effects of global rotation on the evolution of density 
perturbations analytically. 
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3.3.1 Isocurvature perturbations 

Setting DaTZ = in Eq. (32) means that aDaW^ = K^Va and the linear system (30), (31) takes 
the form 



Va = COabV^-{l + w)Za, (33) 
Za = UabZ'' + lKn{l-(^,)Va--^-Da{B''Vh), (34) 



where we have used the background relation = k/x(1 + w)/2 and the barotropic expression 
DaP = {lJ'Cg/a)T>a. The projected divergence of the these equations, multiplied by the back- 
ground scale a, gives the set of equations governing the linear evolution of isocurvature density 
perturbations in a perturbed Godel universe. 

A = -{l + w)Z, (35) 



Z = |kM1-c,')A-^d2A, (36) 

where A = aD"Pa and Z = dD°'Za by definition.^ Note that the scalar A describes local matter 
aggregations and it is the covariant analogue of the standard density contrast 5p/ p. The system 
(35), (36) leads to the following wavelike equation for the evolution of matter aggregations 

A = - |k/x(1 + k;) (1 - c^) A + c^D^ A , (37) 

and subsequently leads to 



\k) = - 



\kii{1 + w){1 - cl) + 



A(fe) , (38) 



assuming the decomposition A = A(^,)(5*''^^ with T)a^{k) = 0. Note that Q'^^'^ arc locally 
defined scalar harmonics with Q^^^ = and Vp'Q^^^ = —{k"^ / a^)Q'^^\ The last equation has an 
oscillatory (i.e. neutrally stable) solution as long as —fi < p < Note that the stability of the 
isocurvature modes is guaranteed even when the fluid pressure vanishes. 

Incorporating DaO;^ in Eq. (34) has increased the overall resistance of the model against 
the gravitational pull of the matter, since gradients in the rotational energy act like pressure 
gradients. Their presence is responsible for the stability of the linear isocurvature perturbations 
demonstrated in (38). 

3.3.2 Perturbations under rigid rotation 

Setting Dau;^ = in Eq. (31) and using thee barotropic fluid relation D^p = {^cl/a)Va, trans- 
forms the system (30), (31) into 

Va = iVabV'' - {1 + w)Za , (39) 
Za = UabZ^-^Kfi{l + 3cl)Va-Yf:^Ba{D^Vb). (40) 



''in deriving Eq. (35) we have used the linear result oD'^Va = A + auiab^'^T^^ ■ An exactly analogous relation 
for the expansion gradients has been used to obtain (36). 
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Multiplying the above with the characteristic background scalar a and then taking their scalar 
parts, as before, we obtain 



A = -{l + w)Z, 
which leads to the following wavelike equation for A(fc) 



A(fe) = 



^k/x(1 + k;)(1 + 3c2)- 



A 



(41) 
(42) 

(43) 



on using the harmonic decomposition A = Aj-^.-jQ^'^) of the previous section. For dust 
(i.e. w = = Cg) there is no pressure support and A grows unimpeded. In the presence of 
pressure, however, there is a effective Jeans length given by 



Aj 



(44) 



since /c//(l + w) = 2co'^ to zero order. Inhomogeneities on scales exceeding Aj collapse under 
the gravitational pull of the matter, but short wavelength perturbations oscillate like sound 
waves. Interestingly, the Jeans scale given above is comparable to Rq, namely to the radius 
of the smallest closed timclike curves (see Eq. (21)). This means that the causal regions of a 
Godel universe containing a fluid with non-zero pressure are stable with respect to linear matter 
aggregations. 

By setting Dqu;^ = we have removed the supporting effect of the centrifugal energy gra- 
dients from Eq. (31). This is the reason for the instability pattern associated with (43). The 
only effect of the Godel background is via the rotational energy contribiition to the overall grav- 
itational pull (see also Sec. 3.2). As a result, the above given Jeans scale is smaller than that 
of the almost-FRW models [53]. Note that the Jeans criterion found here is analogous to that 
obtained in [37]. 



3.4 Directional effects on density perturbations 

The rotation of the Godel universe introduces a preferred direction because of the background 
vorticity vector. The analysis of density perturbations given in the previous sections examines 
the evolution of the scalar A, which describes the average matter aggregation and therefore it 
does not pick out any directional effects. These, however, are incorporated in D^, the projected 
vector that describes variations in the density distribution as seen between two neighbouring 
observers. To reveal the directional effects of the background rotation it helps to implement an 
additional splitting of the Godel space, along and orthogonal to the vorticity vector. As a first 
step, we introduce the unit vector Ua = ojaj^ (where u> = ^/uj^pj^) parallel to uja and use it to 
define the 2-dimensional projection tensor 

fah = Kb - naUb , (45) 

with fab = f(ab), fabn^ = 0, fa'^fb" = fa'' and fa"- = 2. This projects into the 2-D space orthogonal 
to LOa- Employing (45) we define the gradient = fa^^b, with n"Da = 0, and introduce the 
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decomposition = + na'nJ'T^j, orthogonal to and along the rotation axis. Then, the density 
gradients split as 

Va = f>a + VUa , (46) 

where Va = fa'"Db = (a//i)D (j/i is the density perturbation orthogonal to and D — n'''T?a — 
{a/ fi)n^DaiJ. is the density perturbation parallel to the rotation axis. Similarly we write Za = 
Za + Zua, with Za = aD^G and Z = an"-DaO, for the expansion gradients. 

Starting from Eqs. (30) and (31), using the barotropic fluid relation D^p = CgDa/z, expressing 
Dq (D^p) with respect to (DaP) (by changing the order of the covariant derivatives) , and then 
applying the decomposition introduced above, we obtain 

Va = 2(1 - 2cl)eabcio^V'' + i«/x(l + w){2 + c^,)Va + c^B'^Va - 2(1 + w)al)au;^ , (47) 

for density perturbations orthogonal to Ua and 

i> = ^Knil + w){l + 3cl)T> + c^D^P - 2(1 + w)an''BaLo'^ , (48) 

for those along the rotation axis. This reveals a qualitatively different evolution for perturbations 
orthogonal and parallel to uja- The gravitational pull of the matter gets stronger along the 
rotation axis as the fluid sound speed increases above the Cg = 1/2 threshold. The role of the 
background vorticity, on the other hand, is more pronounced orthogonal to the rotation axis 
and depends on the nature of the medium and on the vector product eabc^^V"^- Note that the 
effects of pressure gradients and of gradients in the centrifugal energy density, orthogonal and 
parallel to tOa, are effectively identical. 

One can obtain more quantitative results by looking at certain particular cases. For example, 
consider inhomogeneities parallel to the rotation axis and assume that the vorticity gradients do 
not change along this direction (i.e. set n^DaOJ^ = 0). Then, it is easy to show that V is unstable 
for dust and that there is an associated Jeans length when the fluid has nonzero pressure, 
as predicted in [37]. Alternatively, we may allow for n^DaO;^ 7^ and assume isocurvature 
perturbations (i.e. set aDau'^ = KfiVa). In that case, inhomogeneities parallel to the rotation 
axis are neutrally stable, they oscillate, for all types of matter with w, Cg > 0. Thus, the 
evolution of density perturbations in the direction of the rotation axis is identical to that of 
average matter aggregations. 

The supporting role of the gradient D^o;^ against inhomogeneities orthogonal to the rotation 
axis is also clear, at least for perturbations of the isocurvature type. In this direction the 
background vorticity has an additional effect conveyed by the first term in the right-hand side 
of Eq. (47). This is harder to quantify, however, as it depends in an intricate way on the nature 

of the fluid and on the orientation of the vector eabc'^^V^- Interestingly, the effect of this term 
is reversed as the sound speed of the medium crosses the Cg = 1/2 threshold.^ 

3.5 Vorticity perturbations 

Rotation is controlled by a pair of nonlinear propagation and constraint equations given respec- 
tively by 

'^(a> = -fQ'^a - ^CUrllia + Uab^^ , (49) 

®The projected vector Da contains information about scalar matter aggregations as well as for turbulence and 

distortions in the density distribution of the fluid. One should bo able to isolate and extract this information by 
developing further the decomposition introduced via (45) and (46). This, however, goes beyond the scope of the 
present paper. 
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and 

BaU" = ilaU" . (50) 

The same expressions also hold when we linearise about the Godel background, although then 

only the zero-order vorticity vector contributes to the right-hand side of (49) and (50). For a 
barotropic fluid the vorticity propagation formula takes the linearised form 

Ul^a) =-§(!- fCs) ©t^a + (^ab^^ ■ (51) 

When deriving the above we have used the barotropic expression iia = —[Cg/{iJ, + p)]DaiJ, for the 
4-acceleration, the conservation law (22), and taken into account that on a rotating background 
the projected gradients of scalars do not commute. All these guarantee that curlita = —2c^@uJa- 
Part of the rotational behaviour seen in Eq. (51) is already familiar from studies of the 
perturbed FRW universes. In particular, we see that expansion leads to less rotation when 
Cg < 2/3, but increases it as uja oc a^'^^^ in models with a stiffer equation of state for the 
matter. In the case of contraction the situation is, obviously, reversed - at least until the non- 
linear terms cease to be negligible. The background rotation of the Godel universe, however, 
has introduced additional effects which propagate through the shear term in the right-hand side 
of (51). To get some idea of impact of the shear, we will assume that Ua is a shear eigenvector 
(i.e. set (Tabo;'' = acOa, where a is the associated eigenvalue). Then, Eq. (51) simplifies to 

^{a) = -[1(1- fCs) 0-(T]u;a. (52) 

Accordingly, when a > the extra shear coupling will increase the overall rotation of the model, 
but it will lead to a reduction otherwise. 

To explain the shear effect, recall that positive a means that there is an extra expansion 
along the rotation axis due to the shear effects. Given, the trace-free nature aab, this means 
that there is an overall contraction orthogonal to Ua, which explains why rotation increases when 
a > 0. Clearly, the opposite happens when a is negative. 

3.6 Gravitational-wave perturbations 
3.6.1 Evolution of the shear 

In addition to the vorticity, the shear also affects the linear evolution of gravitational-wave 
perturbations. Here we will only present the relevant propagation and constraint equations, 
which will be used to analyse linear gravity waves within a perturbed Godel universe in the 
following sections. 

For perfect fluid matter, the nonlinear shear evolution is governed by a set of one propagation 
equation and one constraint 

•^(afe) = -|0O-a6 - Eab + D^^jjif,) - a^/^a'^b)'' - OJ(a^^b) + ^(aUb) ! (53) 

D^ab = fDaB + curW„ - 2ea6cW V . (54) 

When linearised about the Godel background, the constraint equation (54) remains the same, 
whereas the propagation equation (53) reduces to 

^ab = -Kb + D(a^f6) - UJ{a(^b) > (55) 
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3.6.2 Evolution of the Weyl components 

A covariant description of the gravitational waves is provided by the electric (Eab) and the 
magnetic (Hab) parts of the Weyl tensor. The two Weyl components support the different 
polarisation states of propagating gravitational radiation and obey evolution and constraint 
equations that are remarkably similar to Maxwell's equations. For a perfect fluid, the nonlinear 
evolution of the two Weyl components is determined by a set of two propagation equations [54] 

E{ab) = -QEab + cur\Hab - ^^(/i + + 2u^ecd(aHb}'^ + ^a^iaEb)'' - u^e^diaEb/ , (56) 

H{ab) = —®Hab — CmlEab + Sac{aHbf — i^'^^cd{aHb}'^ - '^u'^^cdia^b)'^ ^ (57) 

supplemented by the constraints 

D^b = iKD„/x + ea6,a^i^^^-3i^„^.a;^ (58) 
D^Hab = n{n + p)ua-eabcCT\E'''^ + ^Eabu\ (59) 

In addition, the magnetic Weyl tensor is related to the kinematic variables through the nonlinear 
constraint 

Hab = CUrlcTab + B^a^^b) + ^U^a^^b) , (60) 

where curlTa;, = ecd{aB'^T'^b) for symmetric, orthogonally projected tensor Tab- 

Linearised about the unperturbed Godel background, the propagation equations (56) and 
(57) become 

E{ab) = -QEab + CUrlHab - Su^a^^a^b) - (^"^cdiaEb)'^ , (61) 
Hab = -CUllEab - Uj''ecd{aHb)'^ - ^U^^^cia^^b) , (62) 

on using the zero-order relation (4) to express the background Eab tensor with respect to the 
vorticity vector. Similarly, the two constraints (58) and (59) reduce to 

B^Eab = ^KDafi - SHabu'' , (63) 
B'^Hab = K{n + p)LOa + eabc(T''di^''(^'^ + ^EabUj''. (64) 

while (60) does not change. 

3.6.3 Isolating the gravitational-wave perturbations 

Given the symmetric and trace-free nature of Eab and Hab, we can isolate the pure tensor modes, 

namely the gravitational waves, by demanding that Eab and Hab arc also divergence- free at the 
linear level. In practice, this means setting the right-hand side of (63) and (64) to zero. When 
there is no background vorticity, as in perturbed FRW models, we can isolate the gravity waves 
by assuming a barotropic fluid (i.e. setting p = p{fi)) and by switching off both scalar and 
vector perturbations. This is done by introducing the linear constraints D^/i = = DqB = LVa, 
which are consistent to first order [55]- [57]. When the background is rotating, however, the 
aforementioned set of constraints is not enough. Moreover, switching off the model's rotation 
is no longer an option. Here, we are dealing with the rotating Godel background, and we will 
isolate the pure tensor modes by imposing the following self-consistent linear constraints 

BaH = = D„e , D„a;2 = = curlw^ , (65) 

HabUj'^ = = aabUj'^ , Eab'J' = -^^(/^ + P)^a , (66) 
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in addition to the barotropic fluid assumption. Clearly these constraints guarantee the transver- 
sality of Eab and Hab, as well as that of a at, while they maintain a non-zero background rotation. 
Note that constraint (65a) is the standard restriction imposed on perturbed FRW cosmologies, 
while the rest are new. Of the three extra constraints, (65b) guarantees that any perturbations 
in the centrifugal energy that happen to be present are switched off, and that the vorticity vector 
is curl-free. Constraints (66), on the other hand, imply that LOa is an eigenvector of Hab, (^ab and 
Eab to linear order. In the first two cases the corresponding eigenvalues arc zero, whereas in the 
third the eigenvalue is —k{ij, + p)/3. Finally, we point out that the first of the constraints (65a), 
together with the barotropic fluid assumption, guarantees that D^p = 0. This in turn ensures 
that Ua = through the momentum-density conservation. As a result, DaW" = (see Eq. (50)), 
which means that the vorticity vector is solenoidal. 

We check the consistency of any set of constraints by propagating them in time. If every 
constraint is still satisfled, without the need of any extra restrictions, we say that the set is self- 
consistent. Here we can show that the set (65), (66) is self-consistent only when the matter source 
is pressure- free dust (i.e. for p = 0).^ In particular, the consistency of (65a) follows directly 
from the linear propagation equations (30), (31) of the density and the expansion gradients 
respectively. The consistency of the two constraints in (65b) is shown by propagating them in 
time and then by using the linear commutation laws between time derivatives and projected 
covariant derivatives (see [58] for a comprehensive list). Finally, taking the time derivatives of 
HabU^^, CTab^^ and Eaboo^ we can show that constraints (66) are also consistent to linear order. In 
doing so, one needs to employ the zero order relations (4), (10) and use the covariant identities 
given in [58]. 

3.6.4 Linear evolution of gravitational waves 

Once the constraints (65), (66) have been imposed and the pure tensor modes have been isolated, 
Eqs. (61) and (62) reduce to 

Eab = -eEab + CUTlHab-iO%d{aEb)'^, (67) 
Hab = -CmlEab-Lo'^ecdiaHb)'^, (68) 

respectively. In the same situation, the constraint (60) reads 

Hab = cuvlaab + D(aW6) , (69) 

implying that Hab 7^ cuilaab in contrast to the perturbed FRW case. In the Godel spacetime the 
magnetic component of the Weyl tensor vanishes identically (this is a feature that permits the 
existence of the close Newtonian analogue discussed above). Its electric counterpart, however, 
has a nonzero value that is expressed in terms of the model's vorticity (see Eq. (4)). This implies, 
according to [48], that only the gauge-invariance of Hab is guaranteed when linearising about 
a Godel background. Therefore, to avoid any gauge-related ambiguities, we will monitor the 
linear gravity-wave perturbations by looking exclusively into the evolution of Hab- 

To obtain the wave equation for Hab we take the time derivative of Eq. (68) and then employ 
a lengthy, but relatively straightforward, calculation. In the process we apply the constraints 

®In principle, there might be an alternative set of linear constraints, which is less restrictive than (65), (66) 
and still isolates the pure tensor modes of a perturbed Godel universe. It is also possible that a modified set of 
constraints could isolate the linear gravity-wave perturbations, even for matter with nonzero pressure. As yet, 
however, we have not been able to identify any such sets. 
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(65), (66), use the zero order relations (4) and (10), to express E^b and k{ij, + p) with respect to 
the background vorticity vector, and employ the commutation laws and the covariant identities 
of [58]. At the end we arrive at 

Hab-'D^Hab = 0, (70) 

which means that the linear order wave equation of H^b has no vorticity-related source terms, 
despite the fact that the background rotation has not been switched off. Following the above, 
we claim that the Godel universe is neutrally stable against these linearised gravity-wave per- 
turbations. 

4 The case of homogeneous perturbations 

We can also consider the stability of the Godel universe against homogeneous perturbations, 
by ignoring spatial gradients such as Da/x, Da6, DaW^, etc. We will do this by focusing on the 
Raychaudhuri equation, which takes the linearised form 

e = -i(M + 3p)+a;2 + A. (71) 

We can determine the effects of small homogeneous perturbations of the metric that alter the 
matter density and rotation by comparing the relative growth of the vorticity and matter energy- 
density terms in (71). The conservation of angular momentum for rotational perturbations of a 
perfect fluid with equation of state 

p = (7 - 1)^ = 

requires that the angular momentum of an eddy of mass M and comoving scale a be constant; 
thus 

Ma^LO = const. 

Since // oc a~'^'^ and M oc ^a^ we have u> oc a^'^~^ (in agreement with Eq. (52) when a is 
negligible) and so 

— oc a^^-io. 

Thus, for changes that cause expansion (a — ^ 00) the rotation dominates the effects of self gravity 
whenever the fluid state satisflcs 7 > 10/9. This includes the case of radiation (7 = 4/3) but 
excludes that of dust (7 = 1). For perturbations that produce gravitational collapse (a 0) 
the opposite conclusion holds: rotation becomes negligible with respect to the self gravity of the 
fluid whenever 7 > 10/9 but is dominant whenever 7 < 10/9. Superficially, this implies that 
rotation would halt gravitational collapse whenever 7 < 10/9, but the self-gravitating effect 
of the rotation actually contributes to the collapse and a singularity will result so long as the 
matter stresses obey the strong energy condition, in accord with the singularity theorems [4]. 
The detailed behaviour of the rotational collapse depends on the non-linear behaviour of the 
perturbations as a — ^ 0. 

5 Scalar fields in the Godel universe 

Scalar field dominated models are key to contemporary cosmology, given their prominent role in 
inflationary scenarios for the very early universe and as candidates for the dark energy content 
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of the universe today. It is interesting to examine whether the high symmetry of homogeneous 
spacctimes might make them more likely as initial states. If so, they may constraint the existence 
of scalar fields in some way. The symmetries of the Einstein static universe, for example, 
mean that any scalar field that happened to be present will have constant kinetic energy and 
potential [35] . Here we will look into the implications of the Godel symmetries for such a scalar 
field 

Consider a general spacetime filled with a minimally-coupled scalar filed (p and assume that 
Va^i'V"^ < [59]. Then, cp has a perfect fiuid description with respect to a 4- velocity field Ua 
that is normal to the surfaces {4> = constant}. More specifically, the scalar field behaves as a 
perfect fluid with 

fi4, = ^<P^ + V{ct>) and = - V (cj)) , (72) 

relative to Ua = —Va4>/4^ [59]. Note that 4> = — (Va^V"^)^/^ determines the kinetic energy of 
the scalar field and V{(t}) is the associated potential. 

The perfect fiuid description of (j) is necessary if one wants to allow scalar fields in the Godel 
universe. However, since the covariant derivatives of scalars commute, it becomes clear that the 
4-velocity field Ua defined above is irrotational. In other words, the vorticity vector associated 
with Ua vanishes [59], which is in direct contradiction with the generic rotation of the Godel 
spacetime. Putting it in a different way, one can say that the presence of the above Ua implies the 
existence of global spacelike hyper surf aces, something strictly forbidden in the Godel spacetime. 
On these grounds, one could argue that, at least globally, minimally-coupled scalar fields are 
not compatible with the symmetries of the Godel universe. More generally, these considerations 
may offer some insight into the absence of observational evidence for rotation in the universe 
today, to very high precision [60]. If the initial state of the universe is dominated by scalar fields, 
which appear to exist in profusion in string theories [61], then a zero- vorticity initial state will 
be enforced and provide a simple explanation for Mach's 'Principle' even without inflation. Of 
course, any subsequent bout of inflation driven by scalar fields will reduce the effects of pre- 
existing vorticity to levels far below the threshold of detectability today and would be unable 
to generate new vorticity from the scalar fluctuations. As a result, the observation of any large- 
scale vorticity in the universe would be a decisive piece of evidence against an early inflationary 
state [62] and reveal speciflc information about the nature of matter at very high energies. 

It should be emphasised that in this section we have considered an exact Godel spacetime 
containing a single scalar fleld with a timelike gradient (i.e. Va(^V"(^ < 0) and a corresponding 
Segrc type [1, (111)]. However, scalar flelds of different nature are not a priori excluded. For 
example, a scalar fleld with zero gradient (i.e. V = 0); which corresponds to an effective 
cosmological constant and has a Segre type [(1, 111)] is clearly compatible with the symmetries 
of the Godel metric. For more details on the classiflcation of scalar flelds with respect to the 
nature of their gradients the reader is referred to [63]. We should also point out that scalar flelds 
are allowed in rotating Godel- type spacetimes as it has been shown in [44]. 

6 Discussion 

Godel's universe is one of the most intriguing solutions of the Einstein fleld equations. The com- 
mon factor behind its unique and exotic features is the rigid rotation of the model, which is why 
the Godel solution has been widely used to illustrate the possible general-relativistic effects of 
global vorticity and time-travel. In the present article we have looked into the implications of the 
model's generic rotation for the stability of the Godel universe under a variety of perturbations. 
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We first considered scalar-matter aggregations and looked into two different types of inho- 
mogeneities, namely isocurvaturc and perturbations under rigid rotation. We found different 
patterns of stability determined by the presence or absence of gradients in the centrifugal energy. 
These act as effective pressure gradients balancing the gravitational pull of the matter fields. 
Interestingly, the latter is found to have a contribution from the rotational energy as well. This 
is a purely relativistic effect, as opposed to the supporting effect of the vortical energy gradients 
which is Newtonian in nature. When the gradients in the rotational energy are included, as 
occurs for the isocurvaturc perturbations, the model is found to be stable against density inho- 
mogeneities even for pressure-free fluids. In their absence, however, stability is possible only if 
there are pressure gradients and then only on scales below an effective Jeans length. The latter 
is of the order of the largest causal region and smaller than its counterpart in almost-FRW 
cosmologies. This is what happens when dealing with perturbations under rigid rotation. 

The equation of state of the cosmic medium is decisive for the evolution of homogenous 
perturbations. These were studied by comparing the relative growth between the vorticity and 
the matter terms in the Raychaudhuri equation. For perturbations causing expansion, we found 
that rotation dominates over self-gravity when the equation of state obeys p > /i/9, which 
includes the cases of radiation and stiff matter. In the case of contraction the situation is 
reversed. 

The linear rotational behaviour of the perturbed Godel universe was found to follow the 
general pattern familiar from the almost-FRW studies, with some additional effects due to 
shape distortions. We examined these effects by aligning the vorticity vector along the shear 
eigenvectors. We found that, when the shape distortions lead to an extra contraction orthogonal 
to the rotation axis, the vorticity of the model increased, but it decreased otherwise. 

The covariant analysis of linear gravitational waves on rotating backgrounds is complicated 
by the need to select and impose a set a constraints that isolates the pure tensor modes without 
switching off the vorticity. Here, we have introduced for the first time, a set of seven linear 
constraints that does this. We then proceeded to study the gravity-wave evolution by looking 
at the behaviour of the magnetic Weyl tensor. Our choice was based on the fact that this tensor 
has no Newtonian analogue and vanishes in the exact Godel spacetime. The latter property 
guarantees that our analysis is free of any gauge-related ambiguities. The result was a simple 
plane-wave equation for Hab without any vorticity-related source terms, which argues for the 
neutral stability of linear gravity-wave perturbations. 

Finally, we have considered the implications of the Godel symmetries for the presence of 
scalar field. We argue that the absence of any global spacelike hyper surfaces in the Godel 
spacetime, namely the absence of any global irrotational vector field, forbids the introduction of 
scalar fields with timelike gradient vectors. 
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